Effects of hyperons in binary neutron star mergers 
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Numerical simulations for the merger of binary neutron stars are performed in full general rela- 
tivity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and 
neutrino cooling for the first time. It is found that even for the hyperonic EOS, a hypermassive 
neutron star is first formed after the merger for the typical total mass fn 2.7 'Mq, and subsequently 
collapses to a black hole (BH). It is shown that hyperons play a substantial role in the post-merger 
dynamics, torus formation around the BH, and emission of gravitational waves (GWs). In partic- 
ular, the existence of hyperons is imprinted in GWs. Therefore, GW observations will provide a 
potential opportunity to explore the composition of the neutron star matter. 

PACS numbers: 04.25. Dm, 04.30.-w, 04.40.Dg 



Introduction: The properties of the neutron star (NS) 
matter, in particular its equation of state (EOS), are still 
poorly understood. In recent years, the influence of non- 
nucleonic degrees of freedom, such as hyperons, meson 
condensations, and quarks, on properties of the NS mat- 
ter has been discussed extensively [l|, ■ From an ener- 
getic point of view, eventual appearances of these exotic 
phases seem to be natural because of the Pauli princi- 
ple: the chemical potential for neutrons will eventually 
exceed that of exotic particles at some high densities. 
However, it is currently uncertain at which density and 
temperature these exotic particles appear, because their 
interactions are not accurately known. Generally speak- 
ing, presence of such exotic particles results in soften- 
ing of EOS. This effect reduces the maximum mass of 
NS. The recent mass measurement of PSR J1614-2230 
(Mji6i4_223o = 1.97 ± 0.04M Q ) 0] gave a strong impact 
on the properties of the NS matter, indicating that stiff 
EOSs are preferable. However, it does not still constrain 
the existence of exotic particles in NS [2J, |4| . 

There have been a number of numerical studies explor- 
ing effects of such exotic phases on astrophysical phenom- 
ena [oT-[To|. aiming at extracting information of the NS 
matter. In the context of stellar core collapse, effects of 
hadron-quark phase transition on dynamics and neutrino 
signals were studied in 0, Q using the MIT bag model 
for the quark phase. The hadron-quark transition is im- 
printed in the neutrino signals, and in special conditions, 
it could generate a second shock wave which triggers a de- 
layed supernova explosion @. Sumiyoshi et al. Q stud- 
ied the effects of hyperons in failed supernovae adopt- 
ing an EOS developed in [TTj] . By contrast, effects of 
the exotic particles on dynamics and gravitational waves 
(GWs) from binary neutron stars (BNS) mergers and hy- 
permassive neutron stars (HMNS) subsequently formed 
have not been studied in detail (but see recent studies 
in [1, Q), although the coalescence of BNS is one of the 
most promising sources for next-generation kilo-meter- 
size GW detectors [Hj]. 

Among exotic particles, A hyperons are believed to 
appear first in (cold) NS around the rest-mass density 
of p -2-3/w [E EH, where p nuc « 2.8 x 10 14 g/cm 3 
is the nuclear matter density. In this Letter, we 



present the first results of numerical-relativity simula- 
tions for the BNS merger performed incorporating a 
finite-temperature EOS including contributions of A hy- 
perons [1J (Hyp-EOS). We ignore the contributions of 
S hyperons because recent experiments suggest that E 
hyperons feel a repulsive potential [l3[ while A hyperons 
feel an attractive potential [jj) . and hence, A hyperons 
are likely to appear first and to be dominant. In the fol- 
lowing, we report the effects of A hyperons on merger dy- 
namics, black hole (BH) formation, neutrino signals, and 
gravitational waveforms. We show that the emergence of 
A hyperons is imprinted in GWs from the HMNS. 

Setting of numerical simulations: Numerical simula- 
tions in full general relativity are performed using the 
same method and formulation as in Einstein's 
equations are solved in the so-called Baumgarte-Shaprio- 
Shibata-Nakamura (BSSN)-puncture formulation [TtJ 
with dynamical gauge conditions for the lapse function 
and the shift vector; a fourth-order finite differencing in 
space and a fourth-order Runge-Kutta time integration 
are used; a conservative shock capturing scheme with 
third-order accuracy in space and fourth-order accuracy 
in time is employed for solving hydrodynamic equations. 
In addition to the ordinary hydrodynamic equations, we 
solve evolution equations for the neutrino (Y v ), electron 
(Ye), and total lepton (Yi) fractions per baryon, taking 
account of weak interaction processes and neutrino cool- 
ing employing a general relativistic leakage scheme for 
electron neutrinos (v e )-> electron antineutrinos (^ e )j and 
other types (///r) of neutrinos {v x ) [la.[l9j. 

The results with Hyp-EOS are compared with those 
with a finite-temperature nucleonic EOS [20( (Shen- 
EOS). The maximum mass of zero-temperature spherical 
NS for Hyp-EOS is M maXiHyp w 1.8M Q . This is smaller 
than that for Shen-EOS, M max ,shcn ~ 2.2M Q . Although 
Mmax.Hyp is slightly smaller than Mji 6 i 4 _ 2 230, Hyp-EOS 
deserves employing to explore the impact of hyperons on 
the BNS merger, because it can be a viable EOS of the 
NS matter not in extremely high densities and hence for 
studying the evolution of the HMNS that do not have the 
extremely high densities for most of their lifetime. Simu- 
lations with finite-temperature EOSs with exotic phases 
which can produce a stable NS with mass larger than 
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FIG. 1: Maximum rest-mass density, maximum matter tem- 
perature, and maximum mass fraction of hyperons as func- 
tions of time for H135 (solid red), H15 (dotted black), S135 
(dashed green), S15 (short-dotted cyan), and S16 (dashed- 
dotted blue). The vertical thin lines show the time at which 
a BH is formed. 



Mji6i4-2230 should be performed in the future studies. 

Numerical simulations are performed in a non-uniform 
grid [Hdi. The inner domain is composed of a finer 
uniform grid and the outer domain of a coarser nonuni- 
form grid. The grid resolution in the inner zone is chosen 
so that the major diameter of each neutron star in the 
inspiral orbit is covered by 60 and 80 grid points for low- 
and high-resolution runs, respectively, and we confirm 
that the convergence is achieved except for stochastic be- 
haviors due to convective motions. Outer boundaries are 
located in a local wave zone (at ~ 560-600 km along each 
coordinate axis which is longer than gravitational wave- 
lengths in the inspiral phase). During the simulations, we 
check the conservation of baryon rest-mass, total gravita- 
tional mass, and total angular momentum, and find that 
the errors are within 0.2%, 1%, and 2%, respectively, for 
the high-resolution runs. 

This Letter focuses on the merger of equal-mass BNS. 
The gravitational mass of single NS in isolation is Mns = 
1.35 and 1.5M for Hyp-EOS (referred to as H135 and 
HI 5). We perform simulations with the initial condition 
of about 3-4 orbits before the onset of the merger until 
the system relaxes to a quasistationary state. Quasiequi- 
librium states of BNS are prepared as the initial condi- 
tions, as in 21]. We compare the simulation results with 
those for Shen-EOS with M NS = 1.35, 1.5, and 1.6M 
(referred to as S135, S15, and S16) obtained in 0. We 
do not take account of magnetic fields in the present sim- 
ulations, which may play a role in the late-phase evolu- 
tion after the onset of the merger [22| . 

Numerical results: Figure [T] plots the evolution of max- 
imum rest- mass density, /9 max , maximum temperature, 
Tmax, and maximum mass fraction of hyperons, XA, max 
as functions of t — t me rge where i me rge is an approximate 
onset time of the merger. Before the merger (t < t merge ), 
Pmax and T max for H135 and H15 agree well with those 
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FIG. 2: The torus mass as a function of time after the BH 
formation for H135 (solid red), H15 (dotted black), and S16 
(dashed-dotted blue). 
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FIG. 3: Neutrino luminosities for H135 and S135 (top), H15 
and S15 (middle), and S16 (bottom). The dotted vertical 
lines show the time at which a BH is formed. The solid and 
dashed curves correspond to the result with the Hyp-EOS and 
Shen-EOS, respectively. 



for the corresponding Shen-EOS model S135 and S15, be- 
cause XA imax in this phase is small as O(10~ 2 ) and effects 
of hyperons on dynamics are not significant. After the 
merger sets in (i > t meIse ), on the other hand, XA, max 
increases to be > 0.1 in accordance with the increase in 
Pmax, and hyperons play a substantial role in the post- 
merger dynamics and emission of GWs (see below). 

Although the total mass, M, is larger than the max- 
imum mass of the zero-temperature spherical NS for all 
the models, a HMNS is formed after the merger, sup- 
ported by the centrifugal force and thermal contribu- 
tions to the pressure [t| [HI]. They subsequently con- 
tract by emission of GWs, and for H135, H15, and S16, 
they collapse to a black hole (BH) at t — £bh where 
tBH — t 
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11.0, 3.1, and 8.9 ms, respectively. The 
HMNS for S135 and S15 will not collapse to a BH in a 
cooling time, t coo i = E^/L^ ~ 2-3 s, where E t h is total 
thermal energy and L v is total neutrino luminosity [l5l ]. 
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FIG. 4: (a) GWs observed along the axis perpendicular to the orbital plane for the hypothetical distance to the source 
D = 100 Mpc. (b) The effective amplitude of GWs defined by Q.4/|/i(/)| as a function of frequency for D = 100 Mpc. The 
noise amplitudes of a broadband configuration of Advanced Laser Interferometer Gravitational wave Observatories (bLIGO), 
and Large-scale Cryogenic Gravitational wave Telescope (LCGT) are shown together. 
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FIG. 5: /Gw(i) in the HMNS phase, smoothed by a weighted 
spline, for H135 (solid red), S135 (dashed green), and S16 
(dashed-dotted blue). 



After the BH formation for H135, H15, and S16, an ac- 
cretion torus is formed around the BH. Figure [5] plots 
the torus mass, -Mt orus , as a function of time. M tor us in 
a quasistationary state is w 0.082, 0.035, and 0.10A/ Q 
for H135, H15, and S16, respectively, with p max ~ 10 13 
g/cm 3 (see Fig.[T]). For H15, M torus is smaller than those 
for other models because angular momentum is not sig- 
nificantly transported in the outer region of the HMNS 
due to its short lifetime. Interestingly, Mtorus is smaller 
for H135 than that for S16 in spite of the longer lifetime 
of the HMNS. This is likely because the HMNS for HI 35 
just before the BH formation is more compact than that 
for S16. Our results suggest that the emergence of non- 
nucleonic degrees of freedom such as hyperons would have 
a negative impact for driving a short-gamma ray burst 
for which formation of a massive torus is preferred [23| . 

Figure [3] plots neutrino luminosities as functions of 
time. Electron antineutrinos are dominantly emitted for 
all the models as found in the previous studies [HI, [24j . 
For H135, H15, and S16, the /i/t neutrino luminosity 
shows a spike around t ss £bh because high temperature 



is achieved as a result of compression. Although r max is 
very large as > 100 MeV, the peak amplitude is moderate 
because of long diffusion timescale of neutrinos due to the 
high density. Soon after the BH formation, /i/r neutrino 
luminosity steeply decreases because high temperature 
regions are swallowed into the BH, while luminosities of 
electron neutrinos and antineutrinos decrease only grad- 
ually because these neutrinos are emitted via charged- 
current processes from the massive accretion torus [15j . 
These features in neutrino luminosities are quantitatively 
the same for Hyp-EOS and Shen-EOS models, and hence, 
it would be difficult to extract information of the NS mat- 
ter only from the neutrino signal. 

Figure 01a) plots the plus mode (h + ) of GWs as a 
function of i re t — Emerge where t re t is the retarded time, 
extracted from the metric in the local wave zone. The 
GW amplitude is h + < 2 x 10~ 22 for a source at a dis- 
tance D — 100 Mpc. GWs from the inspiral phase (for 
^rct ^ ^mcrgc) agree well with each other for the mod- 
els with Hyp-EOS and Shen-EOS for the same mass. On 
the other hand, quasi-periodic GWs from the HMNS (for 
^rot ^ imerge) show differences. First, the amplitude of 
quasi-periodic GWs damps steeply at the BH formation 
for H135 and H15. This is because the HMNS collapse to 
a BH due to the softening of the EOS before relaxing to a 
stationary spheroid. Second, the characteristic GW fre- 
quency, /gwj increases with time for Hyp-EOS models. 
These facts are clearly observed in the effective amplitude 
[see Fig.gKb)] defined by h eS = 0Af\h(f)\ where h(f) is 
the Fourier transform of h+ — ih x with h x being the 
cross mode and the factor 0.4 comes from taking average 
in terms of direction to the source and rotational axis of 
the HMNS. Reflecting a shorter lifetime of the HMNS 
in Hyp-EOS models, the peak amplitude of h c g(f) is 
smaller, in particular for H15 where the HMNS survives 
only for a short period ~ 3 ms. Reflecting the shift of the 
characteristic frequency, the prominent peak in the GW 
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spectra for Hyp-EOS models (H135 and H15) is broad- 
ened. We here note that the lifetime of the HMNS is 
longer for H135 than that for S16. The reason for this 
broadening of the peak is described as follows in more 
detail. 

In the case that hyperons are absent, the HMNS 
slightly contract simply due to the angular momentum 
loss (weakening centrifugal force) via GW emission. By 
contrast, in the case that hyperons are present, Xa in- 
creases with the contraction of the HMNS, resulting in 
the relative reduction of the pressure. As a result, the 
HMNS further contract. Recent studies showed that /gw 
is associated with the frequenc y of /-mo de which is ap- 
proximately proportional to y/Aln/R^ where Mr and 
i? H are the mass and radius of the HMNS pH [25| . This 
indicates that /gw should increase with time. To see that 
this is indeed the case, we show /gw(= d(f>Np/dt) calcu- 
lated from a Weyl scalar * 4 = * e _ ^ NP in the HMNS 
phase for H135, S135, and S16 in Fig. [5] It is clearly seen 
that the mean value of /gw hi the HMNS phase is ap- 
proximately constant for Shen-EOS models; /gw ~ 2.1 
and 2.5 kHz for S135 and S16, respectively. By contrast, 
/gw for H135 increases with time (from /gw ~ 2.0 kHz 

clt ^ret ^merge — ^ HIS to £^ 2.5 kHz £ttj £ret ^mcrgc — 10 

ms) as the HMNS becomes compact. 

Our results raise a caution that one should be care- 
ful when using the peak frequency of the GW spectrum 
to extract information of the NS matter, because /gw 
evolves with time for Hyp-EOS models. This makes it 
ambiguous to determine /gw an d to relate /gw with the 
HMNS structure. Our results rather suggest a possibil- 
ity that the emergence of hyperons may be captured from 



the evolution of the characteristic frequency of GW and 
the peak width of the GW spectra. Accompanied with 
this broadening, the peak amplitude in h B g decreases by 
a factor of a few. 

Summary: We have reported effects of hyperons on the 
BNS merger in numerical-relativity simulations incorpo- 
rating finite-temperature, both hyperonic and nucleonic 
EOS, microphysical processes, and neutrino cooling. We 
showed that for the adopted hyperonic EOS, a BH is not 
promptly formed and a HMNS is first formed for the typ- 
ical total mass of BNS (~ 2.7M©), as in the nucleonic 
Shen-EOS [3. The HMNS subsequently collapse to a 
BH and a massive torus is formed around the BH. The 
torus mass for the hyperonic EOS is smaller than that for 
the nucleonic EOS. We further showed that for the hy- 
peronic EOS, the characteristic frequency of GWs, /gw> 
from the HMNS increases by a factor of ~ 20-30% during 
their evolution via GW emission, by contrast with that 
for the nucleonic EOS in which /gw is approximately 
constant. Although the hyperonic EOS adopted in this 
Letter has some limitations that it only takes A hyper- 
ons into account and that it cannot produce a stable NS 
with mass Mji6i4-2230i we believe from our results that 
it will be possible to constrain the composition of the NS 
matter via observation of GW from the HMNS. 
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